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Beschreibung 

Verfahren und Vorrichtung zur schnellen Verarbeitung von 
Messdaten mit einer Vielzahl unabhangiger Stichproben 

5 

Die vorliegende Erfindung betrifft ein Verfahren sowie eine 
Vorrichtung zur Verarbeitung von Messdaten, die sich aus meh- 
reren durch zeitlich aufeinander folgende Messungen entstan- 
denen Datensatzen mit einer Vielzahl unabhangiger Stichproben 

0 zusammensetzen, bei dem fur jede im Datensatz enthaltene un- 
abhangige Stichprobe ein gemessener zeitlicher Verlauf, der 
sich aus Messdaten der Stichprobe ergibt, mit dem zeitlichen 
Verlauf zumindest einer in einer Modell -Matrix enthaltenen 
Model lfunkt ion unter Einsatz des General Linear Model vergli- 

5 chen wird, urn das Auftreten bestimmter Muster im Signalver- 
lauf zu uberpriifen. 

In der bevorzugten Anwendung betrifft die Erfindung die Ver- 
arbeitung von Messdaten, die mit einem Messverf ahren der 

0 funktionalen Bildgebung, insbesondere mit der funktionalen 

Magnetresonanztomographie ( f MRI ) , aus einem Obj ektvolumen er- 
fasst werden und sich aus mehreren durch zeitlich aufeinan- 
derfolgende Messungen entstandenen Volumendatensat zen zusam- 
mensetzen, bei dem fur jedes im Ob j ektvolumen erfasste Volu- 

5 menelement ein zeitlicher Signalverlauf , der sich aus Messda- 
ten des Volumenelementes ergibt, mit dem zeitlichen Verlauf 
zumindest einer Model lfunkt ion verglichen wird. 

Insbesondere im Bereich der Medi zintechnik und der medizini- 
0 schen Forschung besteht ein Bedurfnis, Inf ormat ionen liber die 
Gehirnaktivitat in menschlichen und tierischen Organen zu er- 
langen. Die neuronale Aktivierung aufrerst sich in einer Zu- 
nahme des Blutflusses in aktivierten Gehirnarealen, wobei es 
zu einer Abnahme der Blut-Deoxyhamoglobinkonzentrat ion kommt. 
5 Deoxyhamoglobin ist ein paramagnet ischer Stoff, der die Mag- 
net feldhomogenitat verringert und damit mit Hilfe von Magnet- 
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resonanztechniken darstellbar ist, da er die T 2 *-Signal- 
relaxation beschleunigt . 

Eine Legalisation von Gehirnaktivitat wird durch Einsatz ei- 
nes funktionalen Bildgebungsverf ahrens ermoglicht, das die 
Anderung der NMR-Signalrelaxat ion mit einer Zeitver zogerung 
misst . Der biologische Wirkmechanismus ist in der Literatur 
unter dem Namen BOLD-Effekt (Blood Oxygen Level Dependent Ef- 
fect) bekannt . 



Eine schnelle Magnetresonanz-Bildgebung ermoglicht es, den 
BOLD-Effekt in vivo in Abhangigkeit von Aktivierungszustanden 
des Gehirns zu untersuchen. Bei der funktionalen Magnetreso- 
nanztomographie werden hierbei in kurzen zeitlichen Abstanden 
15 Magnetresonanzauf nahmen des zu untersuchenden Ob j ektvolumens , 
beispielsweise des Gehirns eines Patienten, gemacht . Durch 
Vergleich des mit Mitteln der funktionalen Bildgebung gemes- 
senen Signalverlauf s fur jedes Volumenelement des Objektvolu- 
mens mit dem zeitlichen Verlauf einer Modellf unkt ion kann ei- 
20 ne St imulus-spezif ische neuronale Aktivierung detektiert und 
raumlich lokalisiert werden. Ein Stimulus kann z.B. ein soma- 
tosensorischer, akustischer, visueller oder olf aktorischer 
Reiz sowie eine mentale oder motorische Aufgabe sein. Die Mo- 
dellf unktion bzw. die Modellzeitreihe beschreibt die erwarte- 
te Signalanderung des magnetischen Resonanzsignals infolge 
neuronaler Aktivierung. Durch Einsatz schneller Magnetreso- 
nanztechniken, wie beispielsweise der Echo-Planar-Methode, 
konnen geringe zeitliche Abstande zwischen den einzelnen Mes- 
sungen realisiert werden. 

30 

Bei vielen mul t ivariaten stat istischen Analysen wird das so- 
genannte General Linear Model (GLM) fur den Vergleich des ge- 
messenen Signalverlauf s mit dem zeitlichen Verlauf einer Mo- 
dellfunktion eingesetzt. Beim General Linear Model handelt es 
35 sich urn einen Least-Squares-Fit von Messdaten an eine oder 

mehrere Modellf unkt ionen . Mit Hilfe des General Linear Model 
wird ermittelt, welche Linearkombinat ion der Modellf unkt ionen 
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die Messdatenreihen moglichst optimal approximiert . Weiterhin 
kann fur jede Model If unktion berechnet werden, wie signifi- 
kant die Messdaten der Nullhypothese keines Beitrages der je- 
weiligen Model If unktion zur Messdatenreihe widersprechen. Das 

5 General Linear Model wird in vielen Bereichen, wie beispiels- 
weise der Physik oder der Soziologie, zur Analyse von Messda- 
ten eingesetzt. Es eignet sich insbesondere auch zur Analyse 
von Zeitreihen, wie sie bei der funktionalen Magnet resonanz- 
bildgebung ( f MRI ) gemessen werden. Durch Einsatz des General 

0 Linear Model kann hier analysiert werden, ob die gemessenen 
Zeitreihen ein Muster zeigen, das lokaler neuronaler Aktivi- 
tat entspricht. Neben diesem Muster zeigen die Zeitreihen a- 
ber haufig auch andere Charakteristiken, wie beispielsweise 
Drifts oder andere Effekte, die im Rahmen des General Linear 

5 Model ebenfalls "modelliert werden konnen. Dies erlaubt eine 
bessere Analyse der Messdaten als beispielsweise t-Test oder 
Korrelationsverf ahren . So ist es mit dem General Linear Model 
beispielsweise auch moglich, mehrere Effekte im Gehirn paral- 
lel zu analysieren. Auch Gruppenstatist iken uber mehrere Pro- 

0 banden sind moglich. Weitere Anwendungsmoglichkeiten des Ge- 
neral Linear Model finden sich beispielsweise in „Human Brain 
Function" von R. Frackowiak et al., Academic Press. 

Bei der Verarbeitung von Messdaten, die mit dem Verf ahren der 
5 funktionalen Magnetresonanztomographie aus einem Objektvolu- 
men erfasst werden, ist es bisher erf orderlich, die gesamten 
Messdaten, die sich aus mehreren durch zeitlich aufeinander- 
folgende Messungen entstandene Volumendatensat ze zusammenset- 
zen, in den Haupt speicher eines Computers zu laden. Anschlie- 
0 ftend muss fur jedes Volumenelement des vermessenen Objektvo- 
lumens der Signalverlauf bzw. die Zeitreihe aus diesen Mess- 
daten extrahiert und mit der jeweiligen Model If unktion ver- 
glichen werden. So ist bei der bekannten Implement ierung des 
General Linear Model in der frei erhaltlichen SPM-Software 
5 (Wellcome Department of Cognitive Neurology; University of 
London; verof f entlicht unter Gnu Public License; 
http: //www. f il . ion.ucl . ac .uk/spm/) ebenfalls erf orderlich, 
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den kompletten zu analysierenden Datensatz, der bei langen 
f MRI-Studien - eventuell auch uber mehrere Probanden hinweg - 
mehrere 100 Megabytes bis hin zu Gigabytes an Daten umfassen 
kann, in den Haupt speicher des Computers zu laden. Die zu a- 
5 nalysierenden Werte, die zu einer Zeitreihe von Messdaten ge- 
horen, werden extrahiert und das General Linear Model wird 
direkt berechnet . 

Diese bisher realisierte Datenverarbeitung fuhrt daher zu ei- 
Ik 10 nem erheblichen Hauptspeicherbedarf . Da die Messdaten ubli- 
■y^^L cherweise volumenweise , d.h. entsprechend den aufeinanderfol- 

^ B genden Messungen als eine Vielzahl von Volumendatensat zen, im 
Speicher des Rechners vorliegen, fuhrt der Einsatz dieser be- 
kannten Technik auch zu sehr hohen Rechenzeiten, da die ein- 
15 zelnen Zeitreihen liber sehr grofte Bereiche der geladenen Da- 
ten zusammengesammelt werden mussen. 

Eine Vorsortierung der Daten konnte diese Rechenzeit zwar 
vermindern, erfordert jedoch wiederum eine bet rachtliche Re- 
chenzeit und zusatzlichen Speicherbedar f . Zudem kann ein Sor- 
20 tiervorgang erst am Ende der Messung fertig sein, da erst 

dann die Zeitreihen vollstandig vorliegen. Die eigentlichen 
, Berechnungen konnen damit erst am Ende der Messung starten. 

k 

Ausgehend von diesem Stand der Technik besteht die Aufgabe 
■ m5 der vorliegenden Erfindung darin, ein Verfahren und eine Vor- 
richtung zur schnellen und speichersparenden Verarbeitung von 
Messdaten anzugeben, die sich aus mehreren durch zeitlich 
aufeinander folgende Messungen entstandenen Datensatzen mit 
einer Vielzahl unabhangiger Stichproben zusammenset zen . Das 
30 Verfahren soli insbesondere die schnelle und speichersparende 
Berechnung des General Linear Model fur den Fall vieler unab- 
hangiger Stichproben, z.B. Messdatenreihen, wie sie bei der 
fMRI auftreten, ermoglichen. 

35 Die Aufgabe wird mit dem Verfahren und der Vorrichtung gemafi 
den Patentanspruchen 1 und 10 gelost. Vorteilhafte Ausgestal- 
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tungen des Verfahrens sowie der Vorrichtung sind Gegenstand 
der Unteranspruche . 

Das vorliegende Verfahren betrifft allgemein die Verarbeitung 
von Messdaten, die sich aus mehreren durch zeitlich aufeinan- 
der folgende Messungen entstandenen Datensatzen mit einer 
Vielzahl unabhangiger Stichproben zusammensetzen. 

Im Folgenden werden das Verfahren und dessen Vorteile ohne 
Einschrankung auf dieses Anwendungsgebiet anhand der Verar- 
beitung von Messdaten erlautert, die mit einem Messverf ahren 
der funktionalen Bildgebung, insbesondere mit dem Verfahren 
der funktionalen Magnetresonanztomographie , aus einem Objekt- 
volumen erfasst werden und sich aus mehreren durch zeitlich 
auf einander folgende Messungen entstandenen Volumendatensat zen 
zusammensetzen. Bei diesem Verfahren wird fur jedes im Ob- 
jektvolumen erfasste Volumenelement ein zeitlicher Signalver- 
lauf aus den Messdaten dieses Volumenelementes ermittelt und 
mit dem zeitlichen Verlauf zumindest einer in einer Modell- 
Matrix G enthaltenen Model lfunkt ion unter Einsatz des General 
Linear Model verglichen, urn das Auftreten bestimmter Charak- 
teristika oder Muster im Signalverlauf zu uberprufen. Das 
vorliegende Verfahren zeichnet sich vor allem dadurch aus, 
dass fur den Vergleich erf orderliche Berechnungen aus den 
Messdaten in einer sich aus der zeitlichen Abfolge der Mes- 
sungen ergebenden Reihenfolge der Datensatze bzw. Volumenda- 
tensatze jeweils fur alle relevanten Messdaten eines Volumen- 
datensatzes durchgefuhrt und als Zwischenergebnisse gespei- 
chert werden, indem die Zwischenergebnisse aus dem jeweils 
unmittelbar vorangehenden Volumendatensat z mit den neuen Be- 
rechnungen aktualisiert werden, so dass jederzeit Zwischener- 
gebnisse berechnet werden konnen und nach dem einmaligen 
Durchlaufen aller Volumendatensat ze ein Endergebnis vorliegt, 
das eine Aussage uber das Auftreten der Charakteristika oder 
Muster im Signalverlauf zulasst, wobei das fur das Endergeb- 
nis erf orderliche Quadrat eines Fehlervektors aus der Diffe- 
renz des Quadrates eines aus den Messdaten gebildeten Mess- 
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wertvektors und des Quadrates eines aus der Modellf unktion 
berechneten Modellvektors erhalten wird. 

Im Gegensatz zu den bekannten Verfahren des Standes der Tech- 
5 nik werden beim vorliegenden Verfahren somit die Messdaten 
fur den Vergleich nicht nach Zeitreihen bzw. Volumenelemen- 
ten, sondern volumenweise durchlaufen. Jeder Volumendatensat z 
kann nach Durchfuhrung der Berechnungen mit den darin enthal- 
tenen Messdaten und Abspeichern eines Zwischenergebnisses 
0 wieder aus dem Hauptspeicher der Datenverarbeitungsanlage 

verworfen werden. Aus diesem Grunde ist es nicht mehr erf or- 
derlich, die gesamten Messdaten, d.h. alle Volumendatensat ze 
der Messungen, in den Arbeitsspeicher der Datenverarbeitungs- 
anlage zu laden. Es reicht vielmehr aus, lediglich jeweils 
5 die Messdaten eines einzelnen Volumendatensat zes zu laden. 
Aus diesem Grunde kann das vorliegende Verfahren sehr Spei- 
cherplatz sparend, effizient und schnell durchgefuhrt werden, 
da bei der Berechnung jeweils nur auf kleine Datenbereiche 
zugegriffen werden muss. Die Volumendatensat ze mussen hierbei 
0 insgesamt nur einmal volumenweise durchlaufen werden. 

Bei der Losung der obigen Aufgabe wurde insbesondere beruck- 
sichtigt, dass die heutige Computerarchitektur bei der Be- 
rechnung auf grofien Datenbereichen, insbesondere im Gigabyte- 

5 Bereich, nur sehr langsam ist, auf kleinen Datenbereichen je- 
doch Berechnungen mit deutlich hoherer Geschwindigkeit durch- 
fuhren kann. Mit dem vorliegenden Verfahren, bei dem volumen- 
weise immer nur kleine Datenbereiche zur Berechnung einge- 
setzt werden, lasst sich somit eine deutlich hohere Berech- 

0 nungsgeschwindigkeit erzielen. Weiterhin ist der Speicherbe- 
darf, der fur die Verarbeitung der Daten erforderlich ist, 
beim vorliegenden Verfahren nicht mehr von der Anzahl der Vo- 
lumendatensatze, d.h. der Anzahl der Messungen, sondern im 
Wesentlichen nur noch von der Grofie des vermessenen Objektvo- 

5 lumens abhangig. Ein weiterer Vorteil des vorliegenden Ver- 
fahrens besteht darin, dass alle Berechnungen bereits wahrend 
der Messung durchgefuhrt werden konnen. Dadurch werden auch 
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Echtzeitanwendungen sowie die Anzeige eines immer aktuellen 
Zwischenergebnisses moglich. Weiterhin kann die Modellf unkt i - 
on beim vorliegenden Verfahren erst wahrend der Messung defi- 
niert werden und muss nicht a priori bekannt sein. Das vorge- 
schlagene Verfahren ist sehr effizient parallelisierbar , so 
dass mehrere Computer oder Prozessoren zur Beschleunigung der 
Berechnungen eingesetzt werden konnen. 

Das vorliegende Verfahren wird nachfolgend rein beispielhaft 
nochmals anhand eines f MRI-Datensat zes erlautert, lasst sich 
jedoch in gleicher Weise selbstverstandlich auch auf Daten- 
satze anderer funktionaler Bildgebungsver f ahren oder sonsti- 
ger Messdatem mit unabhangigen Stichproben - entsprechend den 
Messdatenreihen der einzelnen Volumenelemente bei der f MRI - 
anwenden . 

Bei der funktionalen Magnetresonanztomographie wird das Ob- 
jektvolumen, beispielsweise das Gehirn eines Probanden, in 
geringen zeitlichen Abstanden mehrmals dreidimensional abge- 
tastet und aus den Rohdaten durch eine Fourier-Transformation 
die gewunschte Bildinf ormation in Voxeln, d.h. Volumenelemen- 
ten, rekonstruiert . Jedem Volumenelement ist hierbei ein ent- 
sprechender Messwert zugeordnet . Die gesamte Messung besteht 
aus mehreren Volumendatensatzen, die die Messdaten der zeit- 
lich kurz hintereinander folgenden Messungen enthalten. Jeder 
Volumendatensatz umfasst die zu einem bestimmten Zeitpunkt 
erfassten Messdaten der einzelnen Voxel des Ob j ektvolumens . 
Fur die Erkennung von Akt ivierungszustanden des vermessenen 
Gehirns werden die aus den zeitlich abfolgenden Messwerten 
jedes Voxels resultierenden Signalverlauf e mit einer oder 
mehreren Modellf unktionen verglichen, urn den Grad einer Uber- 
einstimmung feststellen zu konnen. Gerade beim Einsatz des 
verbreiteten General Linear Model (GLM) gemaft dem vorliegen- 
den Verfahren kann dieser Vergleich sehr Speicherplat z spa- 
rend durchgefuhrt werden. Beim General Linear Model wird eine 
Design- bzw. Modell -Matrix G mit einer oder mehreren Modell- 
f unktionen erstellt, aus der in Verbindung mit den Messdaten 
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jedes Volumenelementes Werte, insbesondere eine t-Statistik, 
berechnet werden konnen, aus denen der Grad einer Uberein- 
stimmung mit den Model lfunktionen der Modell -Matrix ersicht- 
lich ist. 

Bei der Durchfuhrung des Verfahrens unter Einsatz des General 
Linear Model werden die Volumendatensatze der Messungen nur 
einmal durchlaufen. Nach Durchlaufen jedes Datensatzes mit 
den dabei durchgef uhrten Berechnungen wird ein Zwischenergeb- 
nis abgespeichert , das beim Durchlaufen des unmittelbar dar- 
auffolgenden Volumendatensat zes durch die neuen Berechnungen 
- aufgrund der neuen Messdaten - aktualisiert und als Zwi- 
schenergebnis dieses neuen Volumendatensat zes abgespeichert 
wird. Durch diese volumenweise Aktualisierung der Zwischener- 
gebnisse ist es moglich, jederzeit ein Ergebnis zu erhalten, 
aus dem beispielsweise die t-Statistik berechnet werden kann. 
Das hierfur erf orderliche Quadrat des Fehlervektors wird aus 
der Differenz des Quadrates des aus den Messdaten gebildeten 
Messwertvektors und des Quadrates des aus der Model lfunkt ion 
berechneten Modellvektors erhalten. Diese einfache Berechnung 
ist moglich, da der Fehlervektor und der Modellvektor auf- 
grund der Least Squares Bedingung senkrecht aufeinander ste- 
hen und addiert den Messwertvektor ergeben. Fur die Berechung 
des Quadrates des Fehlervektors miissen daher nur die Lange 
des Messwertvektors und die Lange des Modellvektors ermittelt 
werden. Das Quadrat des Fehlervektors kann dann uber den Satz 
von Phythagoras berechnet werden. 

Der Modellvektor enthalt dabei die aus der Modellf unkt ion be- 
rechneten Werte, die den Messwerten moglichst nahe kommen 
sollen . 

Das vorliegende Verfahren beruht somit auf einem Update - 
Vorgang, bei dem die Volumendatensatze nur einmal durchlaufen 
werden. Am Ende des Update -Vorgangs stehen die aus den Mess- 
werten fur die Berechnung der t-Statistik er f orderlichen Da- 
ten vollstandig zur Verfugung. 



5 




30 
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In einer moglichen Ausf uhrungsf orm kann der Update-Vorgang 
bzw. der Durchlauf durch die Volumendatensat ze gleichzeitig 
mit der Messung erfolgen, wobei die jeweilige Berechnung je- 
weils unmittelbar nach Erhalt eines neuen Volumendatensat zes 
durchgefuhrt wird. Der Update-Vorgang ist daher in dieser 
Ausf uhrungsf orm mit dem Ende der Datenauf nahme bereits abge- 
schlossen . 

Bei ausreichender Rechengeschwindigkeit ist auch eine Echt- 
zeitanwendung des vorliegenden Verfahrens moglich, bei der 
wahrend der Messung regelmaftig Zwischenauswertungen vorgenom- 
men werden, indem der Update-Vorgang unterbrochen und die Be- 
rechnung der t-Statistik durchgefuhrt wird. Anschlieliend kann 
der Update-Vorgang fur die im Folgezeitraum gemessenen Volu- 
mendatensat ze fortgefuhrt werden. 

Bei dem vorliegenden Verfahren werden Zwischenspeicher beno- 
tigt, die fur den Einsatz des General Linear Model eine Grofte 
aufweisen miissen, die folgender Gleichung genugt : 

Grofie des Speichers = Anzahl der Volumenelemente x \ 
(2 x Anzahl der Modellf unkt ionen + 2) x 4 Bytes ~ /* 




Die vier Bytes als Speicherplat z fur eine Gleit kommazahl vom 
Float-Datentyp haben sich als ausreichend fur die Durchfiih- 
rung des vorliegenden Verfahrens herausgestellt . Die Anzahl 
der Volumenelemente kann beispielsweise 64 x 64 x 32 Werte 
umfassen. Dies entspricht der Grofie eines ublichen Abtastda- 
tensatzes eines Ob j ektvolumens in der funktionalen Magnetre- 
sonanztomographie . 

Das vorliegende Verfahren kann beispielsweise in einem Compu- 
ter eines Magnetresonanz-Scanners eingesetzt werden, so dass 
eine schnelle Nachverarbeitung von fMRI-Daten moglich ist. 
Durch geeignete Implementierung des Verfahrens kann die Dar- 
stellung der Vergleichsergebnisse auch in Echtzeit erfolgen, 
solange die Datensatze nicht zu viele Volumenelemente enthal- 
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ten und die Anzahl der Modellf unktionen nicht zu hoch gewahlt 
wird. Insbesondere ist die Berechnung in Echtzeit nicht durch 
die Anzahl der Volumendatensat ze limitiert, da die Berech- 
nungszeit unabhangig von der Anzahl der bereits hinzugef ugten 
5 Datensatze konstant ist. 



Das vorliegende Verfahren wird nachfolgend anhand eines Aus- 
f uhrungsbeispiels in Verbindung mit der Figur 1, die bei- 
spielhaft die prinzipielle Vorgehensweise bei der Durchfuh- 

10 rung des vorliegenden Verfahrens zeigt, nochmals naher erlau- 
tert. Hierbei wird zunachst auf den Einsatz des General Line- 
ar Model zur Auswertung eines fMRI -Experiments eingegangen 

9* und anschlieftend das vorliegende Verfahren unter Einsatz des 
General Linear Model an der Auswertung einer einzelnen Zeit- 

15 reihe erlautert. 



Zunachst wird die Modell-Matrix G definiert, die sich aus 
mehreren Spalten mit Modellf unktionen zusammenset zt . Die 
Zeitreihe besteht beim vorliegenden Beispiel lediglich aui 
2 0 zwanzig Messungen. Fur jedes Volumenelement werden dabei " v --v^ 
zwanzig Messwerte x m erhalten, wobei m = 0 ... max- 1 und max 
= 20. Als Modell-Matrix wird eine Matrix mit funf Modellf unk- 
tionen gewahlt, d.h. eine Matrix mit fiinf Spalten und zwanzig 
Zeilen. 

^5 In der ersten Spalte steht die Modellf unkt ion der neuronalen 
Aktivierung, die beispielsweise durch 

G = (11111000110011001100) T 
gebildet werden kann. 

In der zweiten Spalte ist ein konstanter Anteil in den Mess- 
30 daten reprasent iert : G m ,i = 1. 

Die weiteren drei Spalten dienen als Hochpassf ilter , d.h. sie 
modellieren langsame Drifts wahrend der Messung auf Cosinus- 
basis. Diese Modellf unkt ionen konnen beispielsweise zu 

G m , 2 = cos G 3 = cos G m 4 = cos 

vmaxj V J \ max J 
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gewahlt werden. 



5 



10 



15 



20 




30 



Die Modell -Matrix G m# j beinhaltet also die Sammlung an Modell - 
f unktionen, welche an die Messdaten x m angepasst werden sol- 
len. 

Bei der Durchfiihrung des Vergleichs der Messdaten mit diesen 
Modellf unktionen wird eine Least-Squares-Schatzung der Para- 
meter durchgef uhrt . Hierbei wird ein Vektor b definiert, aus 
dem die Ubereinst immung der Messdatenreihe mit einzelnen Mo- 
dellf unktionen der Modell -Matrix ersichtlich ist . Dieser Vek- 
tor b wird in folgender Weise berechnet : 



Der Vektor b enthalt im vorliegenden Beispiel funf Werte ent- 
sprechend der Anzahl der Spalten bzw. Modellf unktionen der 
Modell-Matrix . Der jeweilige Wert represent iert , wie die ent- 
sprechende Modellf unkt ion skaliert werden muss, um die Mess- 
daten zu approximieren . 

Als nachstes wird eine Fehlerbetrachtung durchgef uhrt , um ei- 
ne Aussage liber die Qualitat der Schatzung zu erhalten. Hier- 
zu wird der Wert a berechnet, der sich in folgender Weise er- 
gibt : 



wobei e = G • b - x . 

Schlieftlich wird die so genannte t-Statistik gerechnet, bei 
der fur jedes Volumenelement ein Wert t erhalten wird, dessen 
GroBe ein Mali fur den Grad der Ubereinstimmung der Messdaten- 
reihe mit der betrachteten Modellf unkt ion darstellt . Bei 
funktionalen Bildgebungsverf ahren wie der hier beispielhaft 
erlauterten funktionalen Magnetresonanztomographie ist es ub- 
lich, diese t-Werte zu visualisieren . Fur die Berechnung der 
t-Statistik muss zunachst ein Kontrast definiert werden, der 
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die interessierende Modellf unkt ion angibt . Im vorliegenden 
Fall kann dieser Kontrast c durch folgenden. Wert gegeben 
sein: c = (10000) T . 



Der t-Wert ergibt sich dann in folgender Weise 



t := 



cb 



Bei den bisher im Stand der Technik durchgef uhrten Implemen- 
tierungen wird aus dem f MRI -Datensat z nacheinander die Zeit- 
reihe jedes Voxels extrahiert und der obigen Analyse unterzo- 
gen . Hierfiir ist es jedoch erf orderlich, fur die Berechnung 
die gesamten Messdaten alle Volumendatensat ze gleichzeitig in 
den Arbeitsspeicher eines Rechners zu laden. Neben dem hohen 
Speicherplat zbedarf fiihrt dies auch zu einer entsprechend 
langsamen Berechnungsgeschwindigkeit . 

Demgegeniiber konnen mit dem vorliegenden Verfahren, wie es in 
der Figur 1 in einer Ausf iihrungsf orm dargestellt ist, sowohl 
der Speicherplatzbedarf verringert als auch die Rechenge- 
schwindigkeit erhoht werden. Bei dem vorgeschlagenen Verfah- 
ren wird beim Update -Vorgang der f MRI -Datensat z volumenweise , 
d.h. ein Volumendatensat z nach dem anderen, durchlaufen und 
es werden in jedem Volumendatensat z fur jedes Volumenelement 
die jeweiligen Berechnungen durchgef uhrt . 

Dies kann in folgender Weise erfolgen. Beim Update -Vorgang , 
bei dem alle Volumendatensat ze einmal durchlaufen werden, 
wird fur jeden neuen Volumendatensat z und jeden Voxel der 
Vektor GTx aktualisiert , wobei jeweils nur auf die neue Zeile 
der Modell -Matrix und auf den neuen Messwert zugegriffen 
wird : 

GTxi = GTXi + G m ,i • x m . 

GTx bezeichnet hierbei den Vektor G T • x. Diese Berechnung 
wird in jedem neuen Volumendatensat z fur jedes Volumenelement 
durchgefuhrt und kann auch voxelweise auf verschiedene Compu- 
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ter bzw. Prozessoren verteilt werden. Hierdurch wird der vom 
vorhergehenden Volumendatensatz erhaltene Wert fur die ein- 
zelnen Komponenten des Vektors GTx jeweils aktualisiert . Es 
versteht sich von selbst, dass dieser Wert bei Durchfuhrung 
5 der Berechnung fur den allerersten Volumendatensatz auf Null 
gesetzt werden muss. 

Die Durchfuhrung dieser Berechnung fur jeden Volumendatensatz 
impliziert, dass ein ausreichender Zwischenspeicher fur den 
jeweiligen Zwischenwert von GTx vorhanden sein muss, der der 
10 Anzahl der Volumenelemente im Volumendatensatz mult ipliziert 
mit der Anzahl der Modellf unktionen in der Modell -Matrix ent- 
spricht. 

Weiterhin wird bei diesem Durchlaufen der Datensatze fur je- 
15 den neuen Volumendatensatz und jeden Voxel das Quadrat XX des 
Messwertvektors aktualisiert, wobei jeweils nur auf den neuen 
Messwert zugegriffen wird: 



Gleichzeitig kann bei diesem Durchlauf durch die Volumenda- 
tensatze fur jeden neuen Volumendatensatz auch die Matrix GTG 
(entspricht G T • G ) durch folgende Vorschrift aktualisiert wer- 
den : 



Selbstverstandlich kann diese Berechnung, da sie keine Mess- 
daten enthalt, auch in einem. Schritt vor dem Durchlauf, wie 
in der Figur gestrichelt angedeutet, oder nach diesem Durch- 

30 lauf durchgefuhrt werden. Bei einer Modell-Mat rix, deren Mo- 
dellf unktionen erst im Laufe der Messung definiert oder ver- 
andert werden, oder wenn Zwischenergebnisse berechnet werden 
sollen, ist jedoch die Durchfuhrung gemaft der obigen Vor- 
schrift, d.h. die volumenweise Berechnung wahrend des Durch- 

35 laufes, erf orderlich . 



XX = XX + (x m ) 2 



20 




GTG i/k = GTG ifk + G m ,i • G m , k . 



V 
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Nach Erhalt der Matrix GTG wird diese invertiert. Dies kann, 
da diese Matrix eine kleine symmetrische, reelle Matrix ist, 
im vorliegenden Beispiel eine 5x5 Matrix, beispielsweise 
mittels einer LU-Dekomposition erfolgen. Selbstverstandlich 
5 sind jedoch auch andere Invertierungsalgorithmen einsetzbar. 
Schlieftlich wird fur jedes Volumenelement der Vektor b be- 
rechnet, der sich durch 

b = GTG" 1 • GTx 

10 ergibt. 

Nach Durchfiihrung des Update-Vorgangs liegen somit die Werte 
fur den Vektor b und die Pseudoinverse (G T • G)" 1 f est . 



W 

r ' 

Zur Berechnung des Wertes fur das Quadrat der Fehlervektors 
15 e, im Folgenden als EE bezeichnet, wird ausgenutzt, dass der 
durch das Modell darstellbare Teil des Messwertevektors , als 
Modellvektor M bezeichnet, und der Vektor E aufgrund der Le- 
ast Squares Bedingung senkrecht aufeinander stehen und ad- 
diert den Messwertevektor X ergeben. Somit miissen lediglich 
20 die Lange des Messwertevektors und die Lange des Modellverk- 
tors ermittelt werden, urn daraus uber den Satz von Phythago- 
ras EE berechnen zu konnen. 

Hierbei gilt: 

^ M = G • b MM = M • M 

max-1 

MM = (G • b) • (G • b) = q • q = £q. -q. 

i=0 

sp-\ 

30 mit: qi = Y* G U b j 



Eingesetzt ergibt sich; 



max- 1 sp-\ sp-l 

mm = 2 E^/VE^'^ 

i=0 j=0 k=0 
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sp-\ sp-l 1 

£Efo./G i > j -b k 

7=0 i=0 J 

Die innere Klammer stellt die Matrix GTG" 1 dar, die bereits 
berechnet wurde . Somit gilt: 

MM - MM + GTGj.k-bj-bk 

so dass MM schnell berechnet werden kann und direkt zu EE 
fuhrt: 

EE = XX - MM. 

Fur die bisher durchgef uhrten Berechnungen gilt: 

i = 0 . . sp-l , 
k = 0 . . sp-l , 
j = 0 . . sp-l, 
m = 0 . . max- 1 , 

wobei max der Anzahl der Messungen bzw. Volumendatensat ze und 
sp der Anzahl der Spalten der Modell -Matrix G entsprechen. 

Nach jedem Hinzufiigen eines Volumendatensat zes oder auch nur 
am Ende der Messung kann nun die t-Statistik fur diesen Ver- 
gleich berechnet werden. Hierbei wird im vorliegenden Bei- 
spiel zunachst ein Skalierungsf aktor Scale fur den t-Wert 
nach f olgender Vorschrif t berechnet : 

Scale = Scale + ci ■ GTG^j • Cj # 

wobei c dem Kontrastvektor entspricht. 

Diese Berechnung wird fur jedes Volumenelement durchgef uhrt . 
Weiterhin wird fur jedes Volumenelement das a 2 berechnet, das 
sich nach f olgender bekannter Vorschrif t ergibt : 



max- rang(G) 



i=0 



Aus diesem Wert kann der sogenannte t-Wert durch 

t + ci ■ bi 

fur jeden Voxel getrennt berechnet werden. Dieses t wird 
schlieftlich noch durch 4cr2 Scale dividiert, urn den endgulti- 
gen t-Wert fur jedes Voxel zu erhalten, der beispielsweise 
durch Uberlagerung auf eine herkommliche Magnetresonanzauf- 
nahme des Obj ektvolumens visualisiert werden kann. Die letz- 
ten Berechnungen konnen ebenfalls voxelweise auf unterschied 
liche parallel rechnende Computer bzw. Prozessoren verteilt 
werden, urn die Berechnung zu beschleunigen . 

Figur 1 zeigt weiterhin die durch die gestrichelten Pfeile 
angedeutete Moglichkeit, wahrend des Durchlaufs der Datensat 
ze ein Zwischenergebnis zu berechnen und auszugeben. 
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Patentanspriiche 

1. Verfahren zur Verarbeitung von Messdaten, die sich aus 
mehreren durch zeitlich aufeinander folgende Messungen ent- 
5 standenen Datensatzen mit einer Vielzahl unabhangiger Stich- 
proben zusammenset zen, bei dem fur jede im Datensatz enthal- 
tene unabhangige Stichprobe ein gemessener zeitlicher Ver- 
lauf, der sich aus Messdaten x der Stichprobe ergibt, mit dem 
zeitlichen Verlauf zumindest einer in einer Modell -Matrix G 
10 enthaltenen Model lfunkt ion unter Einsatz des General Linear 
Model verglichen wird, um das Auftreten bestimmter Charakte- 
ristika im gemessenen zeitlichen Verlauf zu uberprufen, 
dadurch gekennzeichnet, 
dass fur den Vergleich er f orderliche Berechnungen aus den 
15 Messdaten in einer sich aus der zeitlichen Abfolge der Mes- 
sungen ergebenden Reihenfolge der Datensatze jeweils fur alle 
relevanten Messdaten eines Datensatzes durchgefiihrt und als 
Zwischenergebnisse gespeichert werden, indem die Zwischener- 
gebnisse des jeweils unmittelbar vorangehenden Datensatzes 
20 mit den neuen Berechnungen aktualisiert werden, so dass je- 
derzeit ein Zwischenergebnis und nach dem einmaligen Durch- 
laufen aller Datensatze ein Ender'gebnis vorliegt, aus dem ei- 
ne Aussage uber das Auftreten der Charakteristika im gemesse- 
nen Verlauf abgeleitet werden kann, 
L 2<5 wobei das fur das Zwischen- oder Endergebnis erf orderliche 
Quadrat eines Fehlervektors aus der Differenz des Quadrats 
eines aus den Messdaten gebildeten Messwertvektors und des 
Quadrats eines aus der Model lfunkt ion berechneten Modellvek- 
tors erhalten wird. 



30 



35 



2 . Verfahren nach Anspruch 1 , 
dadurch gekennzeichnet, 

dass beim Durchlaufen der Datensatze fur jeden Datensatz m 
und jede Stichprobe die Vektorelemente 

GTXi = GTxi alt + G m ,i-G m , k 
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und das Quadrat des Messwertvektors 

XX = XX alt + (x ra ) 2 

5 aus den berechneten Vektorelementen GTxi alt und dem Quadrat 

des Vektors XX alt des unmittelbar vorangegangenen Datensatzes 
berechnet und fur den nachf olgenden Datensatz als Zwischener- 
gebnis gespeichert werden, 

wobei m = 0 . . max- 1 , i = 0 . . sp-1, k = 0 . . sp-1, sp der 
10 Anzahl der Spalten der Modell -Matrix G und max der Gesamtan- 
zahl der Datensatze entspricht, 

und wobei anschlieftend das Quadrat des Modellvektors 



15 



30 



MM = MM alt + GTG j/k -bj-b k 



und der Wert des Quadrats des Fehlervektors 



EE = XX - MM 



2 0 berechnet werden, wobei j = 0 . . sp-1 und bj, b k Vektorelemen- 
te des Vektors b = GTG" 1 ■ GTx sind, die in einem Zwischen- 
schritt berechnet werden. 



3 . Verf ahren nach Anspruch 2 , 
dadurch gekennzeichnet, 

dass beim Durchlaufen der Datensatze weiterhin fur jeden Da- 
tensatz m die Matrixelemente 

GTGi # k = GTGi / j c alt + G m ,i«G m ,k 

aus den berechneten Vektorelementen GTGi, k alt des unmittelbar 
vorangegangenen Datensatzes berechnet und fur den nachfolgen- 
den Datensatz als Zwischenergebnis gespeichert werden. 



35 



4 . Verf ahren nach Anspruch 2 , 
dadurch gekennzeichnet, 

dass vor oder nach dem Durchlaufen der Datensatze die Matrix- 
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elemente GTG i/k in einem Schritt aus der Modell -Matrix G be- 
rechnet werden. 

5. Verfahren nach einem der Anspruche 2 bis 4, 
5 dadurch g e k e n n z e i c h n e t , 

dass nach Beendigung des Durchlaufs aus den berechneten Wer- 
ten fur jede Stichprobe ein t-Wert berechnet wird. 

6. Verfahren nach einem der Anspruche 1 bis 5, 
10 dadurch gekennzeichnet, 

dass fiir die Berechnungen jeweils nur ein Datensatz in den 
Arbeitsspeicher einer Datenverarbeitungsanlage geladen und 
^ nach Durchfiihrung der Berechnungen fur diesen Datensatz wie- 
der verworfen wird. 



15 



20 



7. Verfahren nach einem der Anspruche 1 bis 6, 
dadurch gekennzeichnet, 

dass das Durchlaufen der Datensatze parallel zu den Messungen 
erf olgt . 



8. Verfahren nach Anspruch 7, 
dadurch gekennzeichnet, 

dass das Durchlaufen der Datensatze nach einer vorgebbaren 
Anzahl von Messungen kurz unterbrochen wird, urn ein bisheri- 
ges Ergebnis der Messungen darzustellen, wobei der Durchlauf 
nach dieser Unterbrechung fortgefuhrt wird. 

9. Verfahren nach einem der Anspruche 1 bis 8 zur Verarbei- 
tung von Messdaten, die mit einem Messverf ahren der funktio- 

30 nalen Bildgebung, insbesondere mit der funktionalen Magnetre- 
sonanztomographie, aus einem Obj ektvolumen erfasst werden, 
wobei die durch zeitlich aufeinander folgende Messungen ent- 
standenen Datensatze Volumendatensat ze darstellen und der fur 
eine im Datensatz enthaltene unabhangige Stichprobe gemessene 

35 zeitliche Verlauf dem zeitlichen Signalverlauf eines im Ob- 
j ektvolumen erfassten Volumenelementes entspricht . 
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10. Vorrichtung zur Verarbeitung von Messdaten, die sich aus 
mehreren durch zeitlich aufeinander folgende Messungen ent- 
standenen Datensatzen mit einer Vielzahl unabhangiger Stich- 
proben zusammensetzen, mit zumindest einem Rechenprozessor 
zur Durchfuhrung von Berechnungen, einem Massenspeicher zur 
Speicherung der gesamten Messdaten und einem Arbeitsspeicher 
zur Aufnahme zumindest eines Datensatzes sowie von Zwischen- 
ergebnissen der Berechnungen, 

dadurch gekennzeichnet, 

dass eine Einrichtung zur Durchfuhrung des Verfahrens nach 
einem der vorangehenden Patentanspruche vorgesehen ist. 

11. Vorrichtung nach Anspruch 10 in einer Magnetresonanz- 
Anlage . 
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Zusammenf as sung 

Verfahren und Vorrichtung zur schnellen Verarbeitung von 
Messdaten bei der funktionalen Bildgebung 

5 

Die vorliegende Erfindung betrifft ein Verfahren sowie eine 
Vorrichtung zur Verarbeitung von Messdaten, die sich aus meh- 
reren durch zeitlich aufeinander folgende Messungen entstan- 
denen Datensatzen mit einer Vielzahl unabhangiger Stichproben 

10 zusammensetzen. Fur einen Vergleich des zeitlichen Verlaufs 
jeder erfassten Stichprobe mit dem zeitlichen Verlauf einer 
Modellf unkt ion unter Einsatz des General Linear Model werden 

F die erf orderlichen Berechnungen in einer sich aus der zeitli- 
chen Abfolge der Messungen ergebenden Reihenfolge der Daten- 

15 satze jeweils datensatzweise, d.h. von Datensatz zu Daten- 

satz, durchgefuhrt und als Zwischenergebnis gespeichert, in- 
dem die Zwischenergebnisse des jeweils unmittelbar vorange- 
henden Datensatzes mit den neuen Berechnungen aktualisiert 
werden . 

20 Mit dem vorliegenden Verfahren lasst sich der Vergleich Spei- 
cherplatz sparend, effizient und schnell berechnen. 

FIG 
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Berechne GTG 



I 



Lade ersten Datensatz 
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Berechne GTxj und XX fur alle 
Voxel des Datensatzes 
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Lade weiteren Datensatz 
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Update von GTxj und XX fur 
alle Voxel 
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Zwischenergebnis 
gewunscht oder alle Datensatze 
bereits durchlaufen?. 

Berechne MM 
Berechne EE=XX-MM 
Berechne a und t 



Darstellung von t am Monitor 



